Mixed triplet and singlet pairing in ultracold multicomponent fermion systems with 
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The symmetry properties of the Cooper pairing problem for multi-component ultra-cold dipolar 
molecular systems are investigated. The dipolar anisotropy provides a natural and robust mechanism 
for both triplet and singlet Cooper pairing to first order in the interaction strength. With a purely 
dipolar interaction, the triplet p z -like polar pairing is the most dominant. A short-range attractive 
interaction can enhance the singlet pairing to be nearly degenerate with the triplet pairing. We 
point out that these two pairing channels can mix by developing a relative phase of ±5, thus 
spontaneously breaking time-reversal symmetry. We also suggest the possibility of such mixing of 
triplet and singlet pairing in other systems. 
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The study of ultracold dipolar molecules has recently 
become a research focus of cold atom physics [lj-Q . The 
prominent feature of the dipolar interaction is its d r 2_ 3z 2- 
type anisotropy when the dipolar moments are aligned 
by an external electric field. Considerable progress has 
been made in studying anisotropic condensation of dipo- 
lar bosons Furthermore, dipolar fermionic sys- 
tems provide an exciting opportunity to study exotic 
anisotropic many-body physics of fermions. Experi- 
mentally, a near quantum-degenerate gas of the dipo- 
lar fermion 40 K- 87 Rb has been achieved A number 
of theoretical works have been done for the anisotropic 
Fermi liquid properties of dipolar Fermi gases i in- 
cluding both singlet particle and collective excitations. 

The dipolar interaction also has important effects in 
the Cooper pairing symmetry as studied in Refs [l4l - [l8j . 
In the single component case, the only possible pairing 
channels are of odd parity. Assuming dipole moments 
along the z-axis, the pairing symmetry is mainly of p z 
with slight hybridization with other odd partial wave 
components. Dipolar molecules can have an internal de- 
gree of freedom arising from the hyperfine configurations 
of the constituent atoms. The electric dipolar interac- 
tion is independent of these internal components which 
will be denoted as spin below. The inter-component in- 
teraction opens up the possibility of both spin singlet and 
triplet pairings for the simplest two-component case. It 
would be interesting to study even richer Cooper pairing 
patterns and the competition among them. 

In this article, we show that the dipolar interaction 
favors Cooper pairing in the triplet channel over the sin- 
glet channel. This is an effect directly arising from the 
anisotropy of the dipolar interaction, and it occurs to first 
order in the interaction strength. In contrast, it does not 
appear in the usual condensed matter triplet pairing sys- 
tems such as superfluid 3 He fl9rl2l| : the spin fluctuation 
mechanism based on the strong ferromagnetic tendency 
in 3 He arises from the repulsive part of the 3 He- 3 He in- 
teraction at second order. For a two-component dipo- 
lar fermion system, we find the dominate pairing in the 
spin triplet p 2 -like channel with the purely dipolar in- 



teraction. It can mix with the singlet s + d r 2_ 3z 2 -pairing 
whose pairing strength is tunable through the short range 
non-dipolar s-wave scattering. The mixing occurs with 
a relative phase of ±5 which breaks time-reversal (TR) 
symmetry spontaneously. Pairing in dipolar Fermi gases 
with more than two components is also discussed. 

Samokhin et al. studied non-uniform mixed parity su- 
perfluid states in the presence of dipolar interactions |22j . 
They considered coupling between singlet and triplet 
channels with zero relative phase only. Kabanov [23| also 
considered recently mixture of singlet and triplet pair- 
ing with zero relative phase. While this work was being 
completed, a study of the competition between triplet 
and singlet pairing in dipolar fermionic systems appeared 
that analyzed some of the cases considered here 24 1 . 

We begin with the two-component dipolar fermionic 
system with the electric dipole moments aligned along 
the z-axis. The dipolar interaction reads V^nifi ~ rb) = 
- \^-r 2 \ 3 ^2 (cos flft-fj), where 0ft_f> 2 is the angle be- 
tween (fi — fa) and the electric field E\ d is the electric 
dipole moment. The anisotropy is manifested in the an- 
gular dependence of the V373 with the form of the second 
Legendre polynomial. The Fourier transform of the 3D 
dipolar interaction, V(k) = ^-/^(costffe), only depends 
on the polar angle of k. The Hamiltonian is written as 

h = £[<£)-m]4(a)cc(£) + _L J2 v(k-k') 



PL(k;q)P a p(k';q) 



k,k' ,q 



(i) 



where e(k) = hk 2 / (2m); /1 is the chemical potential; 
Pap{k;q) = c a (—k + q)cp(k + q) is the pairing opera- 
tor; a, (3 refer to j" and J,. Please note that V(k — k') 
depends on the polar angle of the vector k — k' , not the 
relative angle between k and k! . We define a dimension- 
less parameter to describe the interaction strength as the 
ratio between the characteristic interaction energy and 

the Fermi energy: A = E int /E F = § . 

We only consider uniform pairing states at the mean- 
field level, thus set q — in the pairing interaction in 



Eq. [TJ We define P a p(k) — P a p{k\q — 0) which satis- 
fies P a p{k) = —P/3 a (—k). The pairing operators can be 
decomposed into the spin singlet P S i and triplet chan- 
nels P*f z : P sl (k) = ^[P n (Js) - P n $)], P t z ri (k) = 
±[P n (k) + P n {k)], Pf ri {k) = -^[P n (k) - P U (k)], 
and PUk) = -j=[P n (k) + P u (k)], Psl (k) and P£.(k) 

are even and odd functions of fc respectively; P^ ri de- 
scribes the triplet pairing operators whose total spin is 
the cigcnstate of e p • S pa i r with zero eigenvalue. Using 
these operators, the pairing interaction of Eq. [T] with 
q = can be rewritten as 

H V*r = ^7 E { Vtrifck') [ E P lH$) P tH(k') } 

k,k', fj.=x,y,z 

+ V ai $;kr)Pl(ii)P ai (&) }, (2) 

where V tri<si (k;k') = \{V{k~k')TV{k + k')} . V si {k;k') 
is an even function of both arguments k and k' , while 
Vtri(k; k') is an odd function of both. 

The decoupled mean-field Hamiltonian reads 

where 5Z' fc means summation over half of mo- 
mentum space; £(fc) = e(k) — /i; 'l'(fc) = 
(cf(fc), C4.(fc), c|(— fc), c|(— fc)) T . The mean-field gap func- 
tion is defined as A Q/3 (fc) = i £ & V(£ - fc')(l^*/3(£')l)- 
A Q) g can be decomposed into singlet and triplet chan- 
nels as A Q/3 (fc) = A si (k)ia y ap + A tr i,n(k)(ia^a y ) a p. 
The Bogoliubov quasiparticle spectra become E\ t 2(k) = 
ijs,l+ Af >2 (fe) where Af 2 (fe) are the eigenvalues of the 
positive-definite Hermitian matrix A^(fc)A(A;). Its trace 
satisfies \\ (k) + X 2 2 (k) = |A sl (fc)| 2 + £ M |A t „, M (fc)| 2 . The 
free energy becomes 

f = -- y in[2cosh^i-— y 

k,i— 1,2 k,k' ,a=si,(tri,fi) 

x {A^V-^feOAaCfc)}, (4) 

where V^" tri (fc; fc') is the inverse of the interaction matrix 
defined as £ £ fc , V si , tri {%, k')V s ^ tri (k' , k") = 
The gap equations are expressed as 

f d 3 k' - - 
A irijM (fc) = -j —^Vtri(k;k')K{k')A t riAk'), 

A st (k) = - J ±^ Vsi (k;k')K(k')A si (k'), (5) 

where K(k') = taxih[^Ei(k')]/[2Ei(k')]. Eq. [5] formally 
diverges. It can be regularized following the standard 
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FIG. 1: (Color online) The eigenvalues w'iri for the spin triplet 
</> z -channel and w"f d for the singlet s+<i -channel. The latter 
depends on the short-range s-wave interaction Voo-o. 

procedure explained in Ref. (l3.[2o| by replacing the bare 
interaction VtH,si with the renormalized zero energy ver- 
tex functions T tr i,si- At the level of the Born approxima- 
tion, this regularization is equivalent to just introducing 
an energy cutoff of ±u) for £1 in Eq. O where ui is at the 
scale of the Fermi energy. 

To analyze the dominant pairing instability around T c , 
Eq. [5] is linearized. Considering that the strongest pair- 
ing occurs at the Fermi surface and following the stan- 
dard procedure in Ref. [l4j , we define the eigen-gap func- 
tions cj>{ rl si {k) satisfying 

J ^^Vtri, si (k;k')^ si (k') = wj^^ ^k) (6) 

where iVo — ^rfp is the density of state at the Fermi 

surface; w\ ri si are dimensionless eigenvalues; k, k' are at 
the Fermi surface. We neglect the effect of Fermi surface 
distortion on pairing which is a higher order effect in the 
interaction strength A. Then Eq. [3] is linearized into a 
set of decoupled equations 

•:</,.,„ {1 + w triM ln[(2e^)/(7rfc B T)]} = 0. (7) 

The spherical harmonics decomposition of Vt r i,si{k; k') 
reads 

Y-v tri ,si{k;k') = J2 v r H' i n.^(n fc )yi/ m (n & ) > (8) 

;m 

where Vw- m remains diagonal for m but couples partial 
wave channels with V — 1,1 ±2. Vt r i,si on ly have ma- 
trix elements in odd and even partial wave channels, re- 
spectively. Vii\m has the same expressions of the Landau 
parameters of the dipolar Fermi gases given by Fregoso 
et al. [IH except for an overall minus sign due to pair- 
ing and a trivial overall numerical factor difference. For 
I = I' = m = 0, Voo;o = because the average value of 
the dipolar interaction is zero. 
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We diagonalize the matrix Vw-.m to find the dom- 
inant negative eigenvalues which determine the domi- 
nant pairing channels. Two eigenvectors are found with 
eigenvalues much more negative than other channels. 
One lies in the triplet odd parity sector with dominant 
Pz-wave character with slight hybridization with other 
odd parity channels as same as in the single compo- 
nent case [l3|: <^> z (Slfe ) with the most negative eigen- 
value w z ri = — 3.820A, whose eigenvector is ()) z (Q k ) ~ 
0.993Yio — O.I2OY30. The other one lies in the even parity 
spin singlet channel. For the purely dipolar interaction, 
its eigenvalue is w s s f d — —1.935 A and the eigenvector 
lies in the mixed channel of s + d k 2_ 3k 2 as <fi s+d (£l k ) ~ 
O.6Y00 — O.8F20 with nodes. However, this channel is sen- 
sitive to the strength of the short range s-wave scattering, 
which contributes only to the matrix element of Voo;o as 
depicted in Fig. [TJ Experimentally, this scattering can be 
tuned to the scale of the Fermi energy through Feshbach 
resonance, i.e., Vnn:n can become of order 1. As for A, as 
estimated in Ref. [TJ, it could reach 0.1 ~ 0.2. Thus the 
competition between <f> s+d and (j> s can be studied experi- 
mentally in the future. When they become degenerate at 
Voo:o/A ss —3.15, 4> s+d becomes mostly of s-wave charac- 
ter as (j) s+d (n k ) « 0.901Yoo - 0.434F 20 . 

We first consider the case of the dominant triplet pair- 
ing, whose critical temperature is determined from Eq. 
[7]as T Pz ss (2e 7 w/7r) exp(— l/\wf ri \). Its order parameter 
configuration is A tr i^{k) = A^^,* (fifc)^. d is a spin 
space unit vector. Without losing generality, it is taken 
along the z-axis as d = e z . A t „ ; ^(f2 fc ) = A tri e 11 cj> z (flk) 
where 7 is the pairing phase. This phase has line-nodes 
on the equator. It breaks the U c (l) gauge, and the spin 
SU S (2) symmetries, but maintains TR, parity, and a 
^-symmetry of the combined operation of d^ — > — 
and 7 — > 7 + 7r [2|| [27}. Its Goldstone manifold is 
G = [U C {1) <g> SU S (2)}/[U S (1) <g> Z 2 ] = \U c (l) ® S 2 S ]/Z 2 . 
The corresponding low energy excitations include the 
phonon and the spin- wave modes. It supports two differ- 
ent classes of vortices: the usual integer vortex of super- 
fluidity, and the half-integer quantum vortex of superflu- 
idity combined with a 7r-disclination of the d- vector. 

Next we consider the coexistence of the singlet (j> s+d 
pairing and the triplet <p z pairing, when they be- 
come nearly degenerate. The order parameter can be 
chosen as A sl (Q k ) = A Q <f> s+d (n k ) and A triifJ ,(tt k ) = 
A z (j) z (Qk)fifi.z- An important observation is that a rel- 
ative ±-|-phase difference between A S i and A tr i is fa- 
vored, thus spontaneously breaking TR symmetry. This 
can be proved as foll ows: the quasi-particle spectra 
reads E s = VFTAf, with \\ 2 = \A (t> s+d (n k )\ 2 + 
\A z cj) z (tt k )\ 2 ± 2Re(A* z A )(/) s+d (n k )(j) z (fi fe ). The last 
term vanishes for relative phase ±| between An and 
A z . In other words, 4> z + itfi s+d is unitary pairing, i.e., 
A^(fc)A(fc) is an identity matrix up to a factor, and 
A? = A| = lAo^+^fc)! 2 + \A z ^ z (il k )\ 2 \. All other 
phase differences give Xi =/= A 2 , thus are non-unitary 
pairing. To show that unitary pairing is optimal, we 



follow the method presented in Ref. [28| to define the 
function f(x) — —% ln[2 cosh | + x], which satisfies 

£tf{x) > 0, thus /(A?) + f(X 2 ) > 2/(^±% Then the 
first term in Eq. @]is minimized by the unitary pairing, 
and the second term is degenerate for unitary and non- 
unitary pairings. Therefore, <fi z + i<f> s+d is favored. This 
is an exotic fully gapped TR breaking pairing state be- 
cause the nodes of (f> s+d and (j> s do not coincide. It also 
breaks parity but is invariant under the combined parity 
and TR operation. Its Goldstone manifold for the con- 
tinuous symmetry breaking is the same as in the purely 
triplet-^ 2 pairing phase. 

The above analysis can be recaptured in the Ginzburg- 
Landau (GL) framework. The bulk pairing order pa- 
rameters are defined as A^ = ^2 k <f> z (k)A tr i^(k) and 

A s + d = £ fe s + d (fc)A sl (fc). The GL free energy is con- 
structed as 

F = a z {T)Y,\^l\ 2 + a s+d {T)\A s+d \ 2 + P Z \A Z X 
+ /? s+d |A s+rf | 4 + 7iEl A ^ 2 l AS+< T 

+ l2 Y / {^ z ; A 7 AS+dAS+d+cc }^ (9) 

where a z = 7V In (=£-), a s+d = No 111(7^— ), and T s+d 

is defined as T s+d w (2eJu}/n) cxp(-l/|w^ +<i |); (3 Z and 
fis+d terms are the standard quartic terms for the triplet 
4> z and singlet <j) s+d channels, respectively; 71,2 describe 
the coupling between the <f) z and </> s+d -channels. Fol- 
lowing the analysis on a similar problem in Ref. [29j, 
we consider the situation where the two channels are 
nearly degenerate and T Pz is slightly larger than T s + d , 
then the triplet pairing develops first at T c — T p , . Defin- 
ing r = (71 — 2\j 2 \)/(2(3 z ), the condition for the second 
instability to occur is that there exists a lower temper- 
ature V below which \a s+d (T')\ > r\a Pz (T')\. This can 

be satisfied for r < 1, which results in y- = (^ :i ) T=F . 
The ±-|-phase difference between the triplet and singlet 
channel pairing requires that 72 > 0. In Refs. |22l l23j. 
coupling between the singlet and triplet pairings through 
a linear spatial derivative is considered, which leads to 
spatially non-uniform states. Due to spin conservation, 
such a term is not allowed here. 

It is natural to further consider competing pairings for 
even larger number of components represented by the 
internal hyperfine spin degrees of freedom, which is an 
even number. The dipolar interaction is independent 
on them, thus the system has an SU (2N) symmetry. 
The 2N x 2N pairing matrix A a p(k) can be classified 
as N(2N + l)-component symmetric (odd parity) pair- 
ing and N(2N — l)-component antisymmetric (even par- 
ity) pairing, which are generalizations of the triplet and 
singlet channel pairings, respectively. They can be ex- 
plicitly constructed as follows. We define the charge con- 
jugation matrix R as Rij — (— ) l 5i. 2 N-i+i and the time- 
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reversal operators T = RC where C is complex conju- 
gation. For 2N = 2, R reduces to the familiar —ia y . R 
satisfies R T = —R and R 2 = — 1. On the other hand, 
any 2N x 2N Hcrmitian matrix can be expanded in the 
basis of the identity matrix and AN 2 — 1 generators of 
the SU(2N) group. They can be classified as even and 
odd under TR transformation. N(2N + 1) of them are 
TR odd which can be constructed as spin-tensor matri- 
ces P M with odd rank numbers (e.g. spin, spin-octupoles, 
etc.), and N(2N — 1) of them are TR even which can be 
constructed as spin-tensor matrices Q v with even rank 
numbers (e.g. the identity matrix, spin quadrupole, etc.). 
Using TP^T- 1 = -P M , TQ v T-\ = Q v , and R T = -R, 
it can be shown that the matrices of P^R are sym- 
metric and Q V R are antisymmetric, respectively. Thus 
we can decompose A al3 (tt k ) = A aS y jfl (Clk)(Qfj,R)a/3 + 
A svmtV (Q k )(P v R) a 0, where A asy ^(fl k ) (fi = 1 ~ 
N(2N - 1)) and A 3 y m , v (n k ) (u = 1 ~ N(2N + 1)) are 
even and odd functions of fl k , respectively. 

The eigenvalue analysis for competing pairing channels 
is the same as in Eq. [5] by replacing the triplet (singlet) 
pairing with the spin-symmetric (antisymmetric) pair- 
ing, respectively. We next consider the unitary pairing 
in both spin-symmetric and spin asymmetric channels, 
respectively. A convenient choice for the matrix kernels 
is that P z = <j z <S> In for the spin symmetric channel, 
and Q — I2N where In and I2N are the identity matrices 
with N and 2/V-dimensions. The first one corresponds to 
the pairing of J2i=i~2N c l (^) c 2iv-i(~^)> while the second 
corresponds to X«=i~2Ar(~)*~ lc l {k)c\ N _ % {-k) . In the 
^-channel pairing states, the spin-St/ (27V) symmetry is 
broken down into SU(N) x SU(N) x [7(1), thus it has 
2N 2 branches of spin- wave Goldstone modes. The vortex 
configuration is similar to the N = 1 case including the 
usual integer vortex and the half-quantum vortex com- 
bined with a 7r-disclination of spin-texture. Again for the 
mixing between pairing in the (f> s+d (spin antisymmetric) 



and <fi z (spin symmetric) channels, a relative phase ±-| 
is needed to maintain the unitary pairing. 

In summary, we have investigated the competing pair- 
ing symmetries in ultracold multicomponent dipolar 
molecular systems, which provides a wonderful oppor- 
tunity to investigated exotic pairings. We predict that 
the anisotropy of the dipolar interaction provides a well- 
defined pairing mechanism to the spin triplet, or, more 
generally, the spin symmetric channel Cooper pairing. 
The spin singlet even parity channel pairing in the <j> s+d 
channel is tunable by the short range s-wave scattering. 
It mixes with the spin triplet odd parity channel pair- 
ing by developing a relative phase ±| to maintain the 
unitary pairing. This is another type of unconventional 
Cooper pairing breaking TR symmetry. 

We point out that our mechanism of TR breaking mix- 
ing between triplet and singlet pairings is very general. 
For example, in supcrfluid 3 He it was originally proposed 
that pairing would occur in the singlet d-wave chan- 
nel [31(5 induced by the attractive part of the Van der 
Waals interaction. Later, attention focused exclusively 
on pairing in the triplet p-wave channel induced by the 
short-range repulsive interaction [3ll l32j . It is natural 
to expect that both channels could contribute to pair- 
ing at sufficiently low temperatures, leading to coupled 
Balian-Werthamer (BW) triplet pairing and singlet chan- 
nel pairing with a relative phase i^. In metallic super- 
conductors, coupling of an isotropic s-wave state with a 
BW triplet state will lead to a single isotropic gap only 
for the particular case where the relative phase is ±5. 
These possibilities will be discussed separately. 
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0804775, and Sloan Research Foundation. This work was 
also supported in part by NSF under No. PHY05-51164 
at the KITP. 
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